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We investigated a flat multidimensional cosmological model in Gauss-Bonnet gravity in 
presence of a matter in form of perfect fluid. We found analytically new stationary regimes 
(these results are valid for arbitrary number of spatial dimensions) and studied their stability 
by means of numerical recipes in 4 + 1-dimensional case. In the vicinity of the stationary 
regime we discovered numerically another non-singular regime which appears to be periodi- 
cal. Finally, we demonstrated that the presence of matter in form of a perfect fluid lifts some 
constraints on the dynamics of the 4 + 1-dimensional model which have been found earlier. 
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I. INTRODUCTION 



Theories of gravity with the Gauss-Bonnet (GB) term as a correction to the Einstein term 
in the action have been actively investigated for more than 30 years l-jj]. This term appears 
as the first non-Einstein contribution in the extension of the General Relativity (GR) known as 
the Lovelock gravity [(| (the only theory that keeps the equations of motion to be the 2nd order 
differential equations), and string gravity [7]. In 3 + 1 dimensions this term does not contribute to 
the equations of motion, and it can be important either in a combination with another fields (as it 
does in string theory) or in theories containing scalar function of the Gauss-Bonnet invariant (see, 
for example 8]; in this case we, however, have higher-order equations of motion). 

For higher number of spatial dimensions the Gauss-Bonnet term is dynamically important. 
Moreover, in Lovelock gravity (in contrast with string gravity) this term is the only allowed non- 
Einstein contribution in (4 + 1) and (5 + 1) dimensions. This means that it is possible to study 
regimes where this term is not a small corrections to Einstein gravity, but equally important or 
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even dominant. 

In the end of 80th some important results in multidimensional cosmology with the Gauss-Bonnet 
term have been obtained including the analog of Kasner solution in a pure Gauss-Bonnet gravity 
[jj. Further studies reveal interesting differences between (4+l)-dimensional case (this is the lowest 
number of dimensions for GB term to contribute) and higher-dimensional cosmology. This includes 
differences in the form of power-law solutions [jj, [l(3], noncontinuity between vacuum solutions and 
solutions with matter (in (4+1) flat case with non-zero matter density p the limit p — > may 



not coincide with vacuum solutions [12|). Some of the differences can be extended to higher-order 
Lovelock corrections (at least in a flat case [3]). Other example of particularity of (4 + 1) case 
is severe fine-tuning needed for smooth evolution from high density (where GB term dominates) 
to low density (when Einstein gravity is restored) - such evolution requires three of four Hubble 
parameters to be equal with high precision [ijj]. This result have been obtained numerically in the 
theory which includes both GB and Einstein terms, so abovementioned results on dynamics with 
matter [12j (founded for a pure GB gravity) does not directly applicable to this case. Nevertheless, 



as we shell see, nonzero matter density can alter corresponding dynamics significantly. 

The structure of the manuscript is as follows: first, we write down equations of motion for 
the model considered. Then, we analytically investigate stationary case and demonstrate exact 
solutions found. Then we switch to numerical studies and confirm by means of numerical recipes 
the reality of the previously found solutions as well as describe new oscillationary regime. Also we 
demonstrate (also via numerical methods) that the presence of matter in form of a prefect fluid 
lifts some constraints on the dynamics of the model considered. Finally, in the Conclusions we 
summarize the results found. 



II. EQUATIONS OF MOTION 

We consider a flat anisotropic metric in (4+l)-dimensional space-time. We are dealing with 
Einstein-Gauss-Bonnet gravity, and non-vacuum space-time. Lagrangian of this theory have a 
form 



C = R + a£ 2 + Cm, 
where R is Ricci scalar, Cm is the Lagrangian of matter fields and C 2 
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C 2 = R^R^ - tR^BT + R 2 (l) 

is the Gauss-Bonnet Lagrangian. 

We are working in the flat background, so the metric we considering has a form 



gfiu = di a g{-l,a 2 (t),b 2 (t),c 2 (t),d 2 (t)}. (2) 

We use perfect fluid with the equation of state p = wp as a matter source; after varying action 
obtained from (JTJ) using the metric above and perfect fluid as a matter field, one can obtain the 
following equations of motion - there are n dynamical equations and a constraint one. First 
dynamical equation has the form 

2(H b + H 2 ) + 2{H C + H 2 ) + 2(H d + H 2 ) + 2H b H c + 2H b H d + 2H c H d + 

(3) 

+8a 



-in s~i 

(H b + H 2 )H c H d + (H c + H 2 )H b H d + (H d + H 2 )H b H c + ^^ W e {abcd)~^ l+ ^ = 



c 4 

where as usual flj = the rest of them could be obtained via cyclic index permutation. The 
constraint equation has a form 

2H a H b +2H a H c +2H a H d +2H b H c + 2H b H d +2H c H d + 24aH a H b H c H d = -^_ £o (abcd)^ 1+w \ (4) 

III. STATIONARY SOLUTIONS 

The system of equations (J3j) — (OJ) has solutions of a particular type in which all Hubble functions 
are constant (and, correspondingly, scale factors expand or contract exponentially), so, we call 
them as stationary for brevity. Anisotropic exponential solutions in Einstein-Gauss-Bonnet gravity 
with matter were particularly studied in [9] under some assumptions. Here we study maximally 
anisotropic case: 

^ = {-l,e 2 ^,e 2 ^,...,e 2 ^}- 
For solution to be a stationary one it is necessary that right-hand sides of field equations 
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do not depend on time. This means that we need 

-(l + w)J2H i = 0, (5) 

i 

i. e. one of a two conditions should be satisfied: 

1. w = — 1 — matter is a cosmological constant; or 

2. J2i Hi = — volume element is constant. 

As in this section we have some results for an arbitrary number of spatial dimensions, we write 
down the equation of motion for a flat (n + l)-dimensional stationary case (i.e. neglecting all time 
derivatives) : 

8nG 



i<j i<j<k<l 



I - — 

C 4 



(6) 



Oil 



i<k 



- 2 E H i E H k R i + 2H 3 E H i H k+ 



+6J?j ^2 HiH^Hi — 6 ^2 HiHkHiHn 

i<k<l<m 



ijtj k,l^i 
k<l 



8ttG 

—t-WSq. 



i^k 



(7) 



i<k<l 

For w = — 1 we have investigated only (4 + l)-dimensional space (n = 4) with a\ = 1 (a = 02) 
and have obtained the following solutions: 



Hi = H2 = H, Ha = H, 



3 — n A 



AaH' 



where H satisfies condition 



8ttG 1 - 4aH 2 + 16a 2 H 4 
—j~ £ o - — 



l6a 2 H 2 



Hi = His arbitrary, H 2 = H 3 



V-4a' 



1-(Ht/=4^ 8ttG 3 

4 " W^+Aa~H' ~^ £ ° ~ "4a' 



where £ = ±1. 
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rj _rr „ z v ^ 1 v _ -16a 2 H 4 - 1 + WaH 2 z - 4z 

Hl ~ H > H *~W H3 ~-^H> H4 - 4aH(-4aH 2 + 8z + 1) ' 

where 



^ _ 1 - 4aH 2 ± v / 5(4aiJ 2 - l) 2 + 16aF 2 



and i7 satisfies condition 



8vrG 16a 2 tf 4 - 4a# 2 + 1 



c 4 u 16a 2 tf 2 
4. 



H 1 + H 2 (Hi + H 2 ) 2 {8aH l H 2 - 3) + 4ffig 2 
3 ~ 2 y 4(l + 8aF 1 tf 2 ) 

—4aHiH 2 (H 2 + flf ) + i7 3 (Fi + if 2 )(l + 8aH 1 H 2 ) + #i# 2 (l + 4a#i# 2 ) 



(1 + 8aH x H 2 )[H 1 +H 2 + H 3 + l2aHiH 2 Hz) 

and i? 2 satisfies condition 

8vrG _ 4(ff 2 + fl| + HiH 2 )aHiH 2 
~^~ e ° ~ 1 + 8aH 1 H 2 ' 

5. 



H 1 = H 2 = H, H 3 = zH, 

H{288a 2 H 4 z + 8aH 2 z - 2z + 52aH 2 - 1) 

96a 2 H 4 z - 20aH 2 z + 8aF 2 -z-2 + !536zH 6 a 3 + 288a 2 H 4 ' 



/8a# 2 - 2 

where z = — 1 ± W ^ , and H satishes condition 

w 8a# 2 + 1 

8ttG l2H 4 a 



c 4 " £ ° 1 + 8atf 2 ' 
6. 

1 8ttG 3 

H 1 = H 2 = H 3 = ± — iJ 4 is arbitrary, — _ £o = - — . 

V-4a c 4 4a 
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For Y^i Hi = a general solution was investigated for arbitrary number of dimensions. If we 
express 

24 HiHjH k Hi y 11 ■ Y "j Y "'■ Y Hl = 

i<j<k<l i jj^i k^i,j l^£i,j,k 

= Y Hi Y H 3 Y H k(~ H i ~ H 3 ~ H k) = - 3 Y H i Y H i Y Hk = 

= -^Y H "Y H ^-^) = ^Y H rY^ + ^Y H "Y H " = 
= -^ H t+^Y H ^Y H l- H ") = -^Y H t + 3 E H *) 2 

i i j i i 

and similarly for other products of Hubble parameters, we get the equations in the form 

-aicr 2 + 3a 2 (cr| - 2cr 4 ) = 2x , 

( 8 ) 

-aifj 2 + 2a 2 (o"| ~~ 2cJ 4) = (1 + w)*c , 

where 

°"2 = Y H i> a4 = Y H ?> x ° = 

i i 

If Q,\ 7^ and a 2 7^ then we obtain that the general solution is any set of Hi, i = 1, . . . n, that 
satisfies the following conditions: 

sr tt2_ 1 ~ 3w ., 

Oil 



From this we can obtain a particular case of [l| where H\ = H2 = H%, H4 = H§ = . . . = H„ . 
It is interesting to note that from equations ([8]) in the case a\ = and Eq = we have only one 
condition to be satisfied 



and that is in agreement with flo| . It can be also easily seen than the only case for a matter 
other than cosmological constant allowing a stationary solution in the pure Gauss-Bonnet gravity 
is of w = 1/3 type. In this case the condition is 
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One can easily see that (fTU|) turns into ([9]) if £o = 0. These two ([H [TO]) are the only special 
solutions of ([8|) if «i = 0. 



IV. NUMERICAL STUDIES 

Before describing regimes let us make some notes regarding our strategy. First of all, just like 
in previous papers where we numerically studied (4+1)- and (5+l)-dimensional Bianchi-I mod- 



els 
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121 ]. we consider only models with positive initial (and in this section we again use 

dynamical Hi = H~i(t) in contrast with the previous section where due to parametrization used our 
Hubble function were constant). In (4+1) dimensions we specify 3 initial Hubble functions and 
calculate 4th using the constraint equation This expression has a form 

(0) _ ggfgf + 2H<PHV> + 2H^HP - Po 

2ff<°> + +2ffi°> +24aH®HV>HP' 1 J 

and here (and later) we use definition for po (and later for p) to simplify notations: po = WttGeo/c 1 . 
Our goal is to study numerically the influence of non-zero matter density on the evolution scenario. 
We fix 3 Hubble parameters and change p and the equation of state parameter w. In this approach 
two different possibilities exist depending on the sign of denominator. If it is positive, we have no 
restriction on the initial value of matter density p. If the denominator is negative, some maximal 
value of initial density exists (that maximal value satisfy = 0). 

Solving Hi = system with respect to (w, p) gives us values (w s t, p s t) for stationary regime: 



= H 2 a + 4aH b H a H2 + AaH c H b H 2 a + 4aH c H a H<; + H 2 C + H b H a + H c H a + Hj + H c H b 
Wst ~ Hi + HI + HI + H c H b + H c H a + H b H a + l2aH c H b H 2 a + YlaH c H a Hl + 12aH b H a H? ' ( ' 



Pst = -2Hl - 2Hl - 2Hl - 2H c H b - 2H c H a - 2H b H a - 2AaH c H b H 2 a - 24aH c H a Hl - 2AaH b H a H 2 c , (13) 
and that regime corresponds to the one, described analytically in the previous section. 
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Table. Possible trajectories for vacuum (4+l)-dimensional GR+GB model 



From 


To 


Design. 


Avail. 


Standard singularity 


Kasner 


I 


+ 




Recollapse 


II 


+ 




Non-standard sing. 


III 




Non-standard, sing. 


Kasner 


IV 


+ 




Recollapse 


V 






Non-standard sing. 


VI 


+ 



A. Singular vacuum regimes 



We start with listing of possible regimes in (4+l)-dimensional vacuum Bianchi-I model. There 
are two possible outcomes for the past evolution - standard and non-standard singularities, and 
three for the future evolution - non-standard singularity, Kasner regime and recollapse. One can 
naively suppose that all six possible combinations of "past" and "future" regimes are occuring, but 
the reality is much more complicated. We summarized the trajectories in Table. 



As it was noted in type-I regime require fine-tuning of initial Hubble parameters, namely, 
three of them should be equal to each other. Recently [15] it was noted that sometimes the equality 
of two Hubble parameters is enough to achieve this type of regime. 

In contradiction with [14], we failed to find type-Ill transitions. Authors of paper [l^ ] claim they 
found them (around 1.5% of total number of trajectories), so keeping in mind possible numerical 
errors as well as different criteria to stop the integration process we can conclude that this regime 
is very rare if exists at all (see also [ijj for possible reasons behind this). 



B. Regimes in presence of matter 

Since distinguishing of initial conditions with different signs of denominator in Eq. (jlip is 
important for us, we treat these two cases separately. A typical transition diagram on the (p, w) 
plane for the case of positive denominator is given in Fig. [T] (to be precise, it corresponds to the 
vacuum type-II case with positive denominator). White region corresponds to the type-II behavior, 
different grey regions correspond to VI, IV, III and I with increase of darkness of grey; they also 
denoted in figure. 

From this figure we can see that type III evolution (which we have not seen in the vacuum case) 
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I I i I i I i I i I i I 

2 4 6 8 10 

P 



FIG. 1: An example trajectories map in (p,w) coordinates for II case with positive denominator. White 
region corresponds to the type-II transition, different grey regions correspond to VI, IV, III and I with 
increase of darkness of grey; they also denoted in figure. 

and type I evolution (which occurs in the vacuum case only under very specific initial conditions) 
both appear as possible scenarios without any severe fine-tuning of initial conditions. As a result, 
all 6 possible scenarios listed in the Table 1, are possible when matter is taken into account. 

In the transition diagram for the case of negative denominator we have the stationary regime at 
the right limit of the diagram if corresponding w s t lies in the physical zone w s t > — 1. In Fig. [2] as 
a way of example we present a typical example for negative w s t (that is type- VI case with negative 
denominator). The situation at the limit line p = p s t is as follows. 

If w > w s t and w —> w s t + the "lifetime" of the type-II regime is increasing and asymptotically 
reaches infinity when w = w s t- In Fig. [3] we presented the transition in question: Fig. [3ja) repre- 
sents Hubble functions versus time for standard type-II regime (standard singularity to recollapse); 
one can see in Fig. [3^b), where we plot (-Hj(i) — H a (t))/H a (t) that the second singularity (that 
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1 1 I I I I I I I I 

0.2 0.4 0.6 

P 

FIG. 2: An example trajectories map in (p, w) coordinates for VI case with negative denominator. Stationary 
regime occur at w = w s t and p — p s t, i.e. at the intersection of w = w s t with right p-boundary. Stroked 
region in the right-bottom corner denotes periodical trajectories. 

corresponds to recollapse) is isotropic one. With w — > w s t + (and p — > p s t) the situation changes 
according to Fig. 0(c); second singularity is still isotropic (Fig. 03(d)). Finally, when p = p st and 
w = w s t, the stationary regime is reached (Fig. E](e)). 

If w < w s t we have quite different behavior. Namely, it becomes periodical; in Fig. 0] we 
presented it. All through Fig. |3]we used the same value for density p « p s t and decreased equation 
of state w from w ~ w s t — in Fig. Ufa) down to w = — 1 in Fig. Hfe). One can see that 
with decrease of the equation of state the period and amplitude of oscillations are first increasing 
first, and then smoothly turn into quasi-stationary regime at w = —1. This means that small 
homogeneous perturbation of stationary regime does not necessary destroys nonsingular behavior, 
but turn it to periodic oscillations near the stationary point. If p is less enough than p s t then the 
described above periodical regime is not triggered and we have singular regime for any w. 
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FIG. 3: The transition from standard type-II to stationary regime with w st < and w — > w st + 0. Hubble 
functions in the standard type-II regime are presented in (a); relative differences {Hi(t) — H a {t)) / H a {t) in (b) 
demonstrate that the second singularity (that corresponds to recollapse) is isotropic one. With w — > w st + 
the situation changes according to (c): the "lifetime" is increasing with "central" part looks pretty flat; 
the second singularity remains isotropic one (d). Finally, when p = p st and w — w st , the regime becomes 
stationary one (e). 
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FIG. 4: The transition from standard type-II to stationary regime with w st < and w — > w st — 0. We used 
the same p ~ p st for all panels (a)-(e) and decreased the equation of state from w s=s w st — in panel (a) 
down to w = — 1 in panel (e). See text for details. 
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w 




w=w„ 



FIG. 5: An example map of transitions in (p,w) coordinates for 1/3 > w st > case. White corresponds to 
type-II transitions; stroked white periodical trajectories as described in w s t < subsection; light grey — 
type-Ill trajectories; grey - type- VI and dark grey - type-IV. 



In case of w s t > the situation is different from described above. There are no periodic regimes 
in the vicinity of the stationary point. Any small perturbations of initial conditions needed for this 
solution ultimately leads to singular behavior. An example for regime transitions diagram for this 
case is plotted in Fig. [5j In general, we have not seen any periodic oscillations numerically with 
w > 0. 



V. CONCLUSIONS 



We have considered dynamics of a flat (4 + l)-dimensional anisotropic Universe filled with an 
ordinary matter in Gauss-Bonnet gravity. Regarding singular regimes, we have founded that any 
of possible 6 different regimes (with 2 possible initial points in standard or nonstandard singulari- 
ties, and 3 possible future outcomes - recollapse, nonstandard singularity or low-curvature Kasner 
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solution) can be realized without severe fine-tuning of initial conditions. This means that some 
constraints of (4 + l)-dimensional dynamics in Gauss-Bonnet gravity found for vacuum regimes 
are lifter when matter is taken into account. Full description of all transition in both vacuum and 
matter cases is given in 151 ] . 



We have found also two nonsingular regimes - a stationary one with constant values of Hubble 
parameters and zero volume expansion rate (so, ^ Hi = 0) which have been found analytically 
for an arbitrary number of dimensions and oscillatory regime founded numerically in the vicinity 
of the stationary regime for (4 + l)-dimensional case. In our numerical studies we have seen this 
latter regime only for w < 0. 
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